Superluminal tunneling of microwaves in smoothly varying transmission lines 
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Tunneling of microwaves through a smooth barrier in a transmission line is considered. In contrast 
to standard wave barriers, we study the case where the dielectric permittivity is positive, and the 
barrier is caused by the inhomogeneous dielectric profile. It is found that refiectionless, superluminal 
tunneling can take place for waves with a finite spectral width. The consequences of these findings 
are discussed, and an experimental setup testing our predictions is proposed. 
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I. INTRODUCTION 
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£N| | Tunneling is a fundamental process related to the dynamics of various kinds of waves. This phenomenon was already 
pointed out in Gamow's famous work [l| on nuclear a-decay, where the probability of penetration of a-particles with 
' energy E through a potential barrier with height Uq (with Uq > E ) was found to be exponentially small, but finite. 
However, the Gurney and Condon attempt 2\ to find the velocity Vt and the transit time r of such a tunneling revealed 
\^ 1 a basic problem of the theory; how should one define these quantities in the "classically forbidden zone" (Uo > E), 
where both the values of v t and r are imaginary? 

Three decades later the interest in this problem was renewed in Hartmann's paper 0] , where the time of tunneling of 
^ | a particle with energy E through a barrier was determined via the phase of the complex barrier transmission function 
. ■ T = \T\ exp(iip t ), which gives the tunneling time 
c/2 ■ 

^: T = h ^E (1) 

1—1 1 where 2nh is Planck's constant. Using the wellknown expressions for the transmission function of a rectangular 
barrier with width d, Hartmann [3| showed that for the case of a thick barrier (|p| d 3> h, where p is the particle 
momentum), the time r becomes independent of d, and thus for a sufficiently thick barrier the tunneling speed Vt can 
£»o ' reach the superluminal velocity Vt > c. This conclusion, referred to in the literature as the Hartmann paradox, can 
, be deduced from standard formula given in many textbooks. It has stimulated a hot debate which is still going on 
C^l ' 0> S B ■ However, a direct measurement of the electron tunneling time through a quantum barrier has proved to be 
an intricate task. The idea to verify Hartmann's conclusion by means of the classical effect of electromagnetic (EM) 
■ wave tunneling through macroscopic wave barriers was then proposed. Its basis was the formal similarity between 
the stationary Schrodinger equation and the usual wave equation, describing both propagating and evanescent EM 
OO modes in continuous media. 

This idea gave rise to several attempts to examine the tunneling times of EM waves in microwave and optical 
, ranges. Thus, radio wave experiments with tunneling of the TEqi mode in an "undersized" metallic waveguide [H, 
were by some authors interpreted in favor of the concept of superluminal phase times for tunneling EM waves. These 
rS \ conclusions were reached by means of the expression [9j 




tb = \ hrM + ( 2 ) 



generalizing Eq. ([I]) to wave tunneling through a barrier characterized by a complex transmission coefficient. However, 
measurements of the group delay time for the fundamental TEqi mode near the cut-off frequency had shown that 
each of the quantities r (Q} and tb © correspond to experiments only at some restricted frequency ranges. Tunneling 
through such homogeneous opaque barriers is accompanied by strong reflections, attenuation and reshaping of the 
transmitted signal. The destructive interference between the incident and reflected parts of the pulse was considered 
in Ref. (Io| as a mechanism for suppressing the tail of the transmitted pulse, causing a superluminal motion of the 
pulse peak. 

Another approach to this problem is based on the use of a hetero gen eous photonic barrier, characterized by a 
curvilinear profile of the dielectric susceptibility e{z) across the barrier Propagation of waves in this geometry is 
characterized by the following inhomogeneity induced phenomena: 
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a) The appearance of an easily controlled cut-off frequency Q, depending on the gradient and curvature of the profile 

s(z). 

b) Formation of a tunneling regime for waves with frequency to < Q in dielectric media with e > and de/dz < 0. 

c) Reflectionless (non-attenuative) tunneling providing 100% transfer of wave energy through a barrier for some 

frequency lj c < Q. 

These properties were examined in Refs. [ill ll2T | for thin dielectric nanofilms, characterized by a smooth spatial 
variation of e(z) at the subwavelength nanoscale. Access to such films, being available due to recent progress in 
nanotechnology [l3j |. remains, however, a challenging technological problem. 

The purpose of the present paper is to unite the advantages of both microwave and optical tunneling systems, 
to enable measurements of substantial phase shifts of GHz microwaves in the regime of non-attenuative tunneling 
through heterogeneous dielectric layerss, characterized, unlike the optical system, by cm-scales for the inhomogeneity. 
A candidate for such a device is the coaxial transmission line (TL), using the fundamental TEM mode. This TL is 
chosen due to the following salient features, resembling the properties of EM waves in free space [l4[ : 

a) The velocity of the TEM mode is known to coincide with the free space light velocity c. 

b) The polarization structure of the TEM mode contains only the transverse field components of E and H, and not 

any longitudinal component. 

c) The TEM mode does not possess any frequency dispersion. 

The organization of the present paper is as follows: Tunneling of the TEM mode through a segment of a coaxial 
TL, containing a diaphragm made from a gradient dielectric metamaterial is considered in section 2. Conditions 
for reflectionless tunneling of a TEM wave through this diaphragm, accompanied by a substantial phase shift, are 
examined in section 3. The parameters of an interferometric scheme for measurements of this phase shift are presented 
in section 4. Finally, a brief discussion of the results is given in the Conclusion (Section 5). 



II. AN EXACTLY SOLVABLE MODEL FOR AN INHOMOGENEOUS TRANSMISSION LINE 



Let us consider a coaxial waveguide permitting EM wave propagation between two infinitely conducting cylinders 
with radii b and a (b > a) shown in Fig 1. The fundamental TEM mode is assumed to propagate along this 
waveguide (z-direction). The regions z < (region 0) and z > d (region 2) are assumed to be vacuum, whereas region 
1 (0 < z < d) is filled with an inhomogeneous dielectric layer, fabricated from a metamaterial, where the dielectric 
susceptibility is varying such that 

ei^e m W 2 (z), (3) 

where W 2 (z) (with W(0) = 1 and W(d) = 1) is a dimensionless positive function and e m is the maximum value of 
e reached at the ends z = and z — d of the inhomogeneous region. The spatial variations of the current / and the 
voltage V along this transmission line are governed by the wellknown equations [14j 

dl dV dV dl 

+ C^- = and + M^- = (4) 

oz ot oz at 

Here M and C are the inductance and capacitance per unit length. Losses are ignored. The values C n and M n are 
for each region are 

C„ = ^ andM„=^ln^V (5) 



1<) 2^ 

where n — 0, 1, 2 and where £o and /io are the dielectric susceptibility and magnetic permeability of vacuum, respec- 
tively. To solve the system @ in the segment n = 1, we introduce the generating function ip defined by 



Next, considering a harmonic dependence oc exp(— iut) we obtain 



92 nm^= ( 7) 
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where v 2 = l/M n C n such that vq = v 2 = c and Wq = W% — 1. Let us now consider a model where W 2 (z) contains 
two free parameters L\ and L 2 which can be considered as the characteristic spatial scales of the inhomogeneity (see 
Fig 2) 



Wf = 



1 , z z z 



(8) 



These parameters can be expressed in terms of the minimum value of the dielectric susceptibility, £ m i n , reached in 
region 1 due to the profile ©. Introducing the ratio y = L2/2L1 we find [l5| 



£m 1 t d at d 



(9) 



Next we introduce the phase path length rj(z) = J* W(zi)dzi, and the new function ip — F/^/W . This leads to the 
resulting equation 



2 F 
dr] 2 



q 2 F = 



(10) 



with q 2 — (uj 2 / v 2 )(u 2 — 1) and u = Cl/co where u > 1. Here Q is a cut-off frequency arising due to the inhomogeneous 
profile W 2 (z). It is 



n = 2 ^ (1 + ^ 2)1/2 where t* 



(11) 



Unlike the evanescent modes in Lorentz media [17], characterized by a local frequency dispersion of natural media, 
Eq. (11) shows the possibility of tunneling of LF waves (oj < £1 ) in a metamaterial with nonlocal dispersion. The 
solution for F is a sum of a forward and a backward propagating wave, i.e. F = e~ + Qe + such that 



A(e- +Qe+) 



II" 



(12) 



where e = exp(— qrj) and e + = exp(qr)) and A is a normalization constant. Using Eqs. (JSJ) the voltage and current 
are 



Vi = A 



e- + Qe H 



IF 



and h = AqVW 



1 

qL 2 



y 



Li 



}e+) - (e- - Qe^ 



respectively. The continuity of the current and voltage at z = gives us the reflection coefficient 

1 + X 



Rv = 



i-x 



with 



X 



iZ N 
Zi 



1 



l-Q 



2qL 1 1 



where N — \/u 2 — 1 whereas Zq and Z\ are the impedances of regions 0, 2 and 1, respectively, i.e. 

Zq = Z2 

The boundary condition at z = d gives 



^ and Z x = — ^ 



where 



Z* 



= exp(2 gT? o) [pN + l+i] 
[PN-i-i] 

N/3 2f3uy 



Z\ ujLi ^Ji + y 2 

ryo = Vo(d) = L 2 /(l + y 2 ) 1/2 J = Hv+/V-) and V± = i 1 + V 2 ) 1 ' 2 ± V 



(13) 



(14) 



(15) 



(16) 



(17) 



(18) 
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Thus, we have found a spatial structure of the E and H components of the evanescent TEM mode. Unlike the 
traditional concept of a homogeneous coaxial TL, which is known to have no cut-off frequency for the fundamental 
TEM mode [bH ]. our inhomogeneous region for a coaxial waveguide has been shown to provide a cut-off frequency 
for this mode without any deformation of the coaxial cross-section. The drastic consequences of this inhomogeneity 
induced effect for reflectance and transmission of the TEM mode are discussed below. 



III. WINDOWS OF TRANSPARENCY FOR THE EVANESCENT WAVE (REFLECTIONLESS 

TUNNELING OF THE TEM MODE) 

In order to find the reflectance of the gradient wave barrier, we substitute the formula (fT7|) for Q, into Eqs. (fT4|) -(fT5 |) . 
This yields the explicit expression for the complex reflection coefficient 

tanh( Wo ) (l + £■ + (3 2 N 2 ) - 1 (3N 
Rv = 7 ; ^ ■ (!9) 



tan 



HQVo) (l - x - P 2 ^) + lP N + i( W N ~ 7 tanh^o)) ' 



Substitution of (|T9|) into the continuity condition for the voltage at z = 0: Vi n {\ + R v ) — iu>~ A(l + Q), where Vi n is 
the voltage of the incident wave, determines the normalization constant A. Then, substitution of the constant A into 
(|13[) gives the complex transmission function T — \T\ exp(i(pt), i.e. 



with 



|T| = ^l-|i?„| 2 (20) 

f tanh( OT o) (l - £ - P 2 N 2 ) + j/3N \ 
^ t=ta11 \ 2^V- 7 tanh( g , ) j (21) 

The formula (fT9|) gives the reflectance of a single inhomogeneous layer. If the wave is tunneling through m (m > 1) 
neighboring and identical layers, the total reflection coefficient R v and phase (ft are found from and (|2"T]) by the 
replacement 



tanh(q?7o) — > tanh(mgryo) (22) 
We note that the reflection coefficient R v can be zero, which, for a series of m layers, corresponds to the condition 

tanh(mgr? ) = 7 ^v 7 ^" (23) 

(l + £ + e m N 2 ) 

The parameters 7 and 770 arc defined by (|18p . Eq. (|23[) determines the normalized critical frequency u c — Q/lo c that 
results in reflectionless tunneling. 

To optimize the parameter values for an inhomogeneous reflectionless barrier, one has to choose the maximum and 
minimum values of the dielectric susceptibility, e m and £ m i n respectively, and the number of barriers m. The next 
steps are: 

1. Calculation of the values y from ([9]) and vi from (jlip . 

2. Solution of Eq. (|23|) for the unknown normalized critical frequency u c = £1/lu c corresponding to the above values 
for y and v±, 

3. Determination of the product u> c d from the condition following from (jTTJ) 



Ucd= ^+y^ =K (24) 

U c 

As a numerical example of such reflectionless tunneling one can consider the values y% — 1/3, e m — 4.912, m = 2; 
in this case the solution of Eq. is u c — 1.0205 and the constant K = 1.761010 cm/s. Then for a frequency 
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w c /2ir — 1GHz (A = 30cm) we obtain d = 2.8cm <C A. We note that the condition for refiectionless tunneling 
(complete transmission = 0, \T\ = 1) can be fulfilled for different frequencies uj and lengths d, linked by Eq. (|23p . 
An example of the transmission spectrum for refiectionless tunneling for the normalized frequencies u is presented 
in Fig 3a. The transmission is almost complete, \T\ 2 > 0.95 in a finite spectral range (1.02 < u < f.25), forming a 
transparency window for the tunneling mode. 



IV. PHASE EFFECTS IN THE TRANSMITTED MODE: SUPERLUMINAL TUNNELING? 

The phase shift ipt of the tunneling mode ([2Tj) is depicted in Fig 3b, as a function of a normalized frequency u, for 
fixed parameters e m and e m in- For the refiectionless tunneling under discussion this shift is <pt = 1-6 rad. This graph 
can be used for different frequencies u> and lengths d, linked by a constant product Lud. We notice that this constant 
product tud determines the phase shift tpo = ujd/c of a wave with frequency uj accumulated along the path d in free 
space. 

To measure the phase shift <p t of the tunneling wave, and comparing it with (p — ujd/c, one can use the 
interferometer-like scheme depicted in Fig. 4. The TEM mode, produced by a generator G, is splitted between 
two similar coaxial transmission lines. The space between the cylinders in TL 1 is empty, whereas the second one 
contains the inhomogeneous barrier that has been discussed above. After passage through these arms, the two waves 
are interfering, and the relative phase shift Aip = ip t — ifo can be measured. Let us write 

<p t = ( 1 + —) (25) 

and introduce the phase-time delay T and the relevant velocity Vt, linked by the condition v{T = D = t$c. Then, 
putting ipo = ujto and tp t = ujtoVt/c, we find 

— = H , T = and t = - 26 

Thus in the case Atp > (Aip < 0) the scheme discussed will result either in a subluminal (v t < c, T > to) or 
a superluminal (v t > c, T < t ) regime. For the figures related to the above mentioned refiectionless tunneling 
(to = 2, D = 2d) one can find from (|2f p . Aip — 0.43rad, which indicates a superluminal propagation of the evanescent 
mode with vt = f.35c and T = 0.74io- Unlike the phase velocity of propagating wave = w/k, characterizing the 
continuous accumulation of phase (with the wavenumber k real), we consider the velocity connected with the phase of 
evanescent wave, which is not associated with its wavenumber. Some salient features of such a superluminal tunneling 
are: 

f. Unlike the case with the Hartmann geometry [3(, no thick barriers and phase saturation are needed in the 
present case. 

2. The power flow of the TEM mode is completely transmitted by means of the evanescent wave. 

3. The energy density of the evanescent mode inside the barrier (e \E\ 2 + /i \H\ 2 )/2 remains positive in each cross- 
section of the barrier. 

It is interesting to evaluate the group time delay for this type of tunneling by means of equations ([T]) and ([2]). 
Since we consider the case of refiectionless tunneling (|T| = 1), the second term in @ vanishes, and we can use |T]). 
Presenting the derivative d(p t / doj from (|2I|) we can write 

d , , > K , , 

t = — and w(u) = g— 27 

The constant K is given in (|24p . Thus, the time r is proportional to the length d when the normalized frequency u is 
given. For the example discussed, i.e. uj = 1GHz and d = 2.8cm, the group delay is negative, or r = —0.6ns. Here |r| 
is subluminal, i.e. |t| > D/c = 0.185ns. By considering the frequency uj = 4GHz, (d = 0.70cm) keeping the product 
ujd constant one finds r = — 0I5ns, which implies a superluminal value |r| < D/c. 

One has to emphasize that these tunneling wave phenomena do not violate the Einsteinian causality related to 
travelling waves. 



6 



V. CONCLUSIONS 

In conclusion, we stress the following points: 

1. The possibility of reflectionless tunneling of microwave TEM mode in a coaxial waveguide with a gradient 
profile (i.e. gradient media) given by e(z) has here been examined in the framework of an exactly solvable 
model. This model describes the effects of heterogeneity-induced dispersion and the controlled formation of a 
cut-off frequency f2 in the wanted spectral range. Unlike evanescent waves in Lorentz media with a natural local 
dispersion 5], the tunneling of LF modes (lu < Q) arises in gradient media due to non-local dispersion. 

2. In contrast to the traditional concept of tunneling in media with e < 0, another mechanism of tunneling, 
with e > but Ve < 0, has been considered here. The possibility of a superluminal phase shift, arising 
in subwavelength wave barriers (beyond Hartmann's condition \pd 3> h\), is demonstrated. An experimental 
setup, illustrating such tunneling, is suggested. 

3. The results, obtained here for microwaves, can easily be generalized to other types of EM waves, keeping the 
condition of a constant product tod, with the other parameters of tunneling media being the same. Moreover, 
the tunneling effects discussed here seem to be rather general, occurring for different types of waves satisfying 
heterogeneous wave equations [1 71 ] for media with continuous spatial variations of its parameters. 
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Figure Captions 

Fig 1: A schamtic picture of a coaxial waveguide with a gradient wave barrier in segment 1 between the empty 
segments and 2. A spatial distribution of the normalized impedance W 2 (z) plotted as a function of z (1/2 = 1/3) 
can also be seen. 



Fig 2: The transmittance spectra for the TEM mode tunneling through two wave barriers in the coaxial waveguide 
(e m = 4.912, e min = 3.69, m = 2, y 2 = 1/3); 

a) The spectral window of transparency in the range 1 < u < 10 and 1.02 < u < 1.25. 

b) The phase spectrum of the tunneling mode in the range 1 < u < 10 and 1.02 < u < 1.25. 

Fig. 3: An interferometer like scheme for measurement of phase shift of the tunneling mode. G is the generator of 
the TEM mode, 1 is an arm containing two adjoining gradient impedance layers, 2 is an arm containing the empty 
coaxial, P is a phase detecting element. 
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